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Abstract. We address nonlocality of continuous variable systems in the presence 
of dissipation and noisc. Three nonlocality tests have been considered, based 
on the measurement of displaced-parity, field-quadrature and pseudospin-operator, 
respectively. Nonlocality of twin-beam has been investigated, as well as that of its non- 
Gaussian counterparts obtained by inconclusive subtraction of photons. Our results 
indicate that: i) nonlocality of twin-beam is degraded but not destroyed by noise; 
■^j- , ii) photon subtraction enhances nonlocality in the presence of noise, especially in the 

low-energy regime. 
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Qh, 1. Introduction 

+-> 

Nonlocality, i.e. the existence of correlations which cannot be explained by any local 
hidden variable model, is perhaps the most debated implication of quantum mechanics. 
During the last decade other aspects of nonlocality, in addition to generating nonlocal 
correlations, have been discovered. For example, the possibility of teleporting and 
effectively encoding information, as well as the ability to perform certain computations 



exponentially faster than any classical device. 

Realistic implementations of quantum information protocols require the investiga- 
tion of nonlocality properties of quantum states in a noisy environment. In particular, 
the robustness of nonlocality should be addressed, as well as the design of protocols to 
preserve and possibly enhance nonlocality in the presence of noise. 

The evolution of nonlocality for a twin-beam state of radiation (TWB) in a thermal 
environment was studied in Ref. [T| by means of the displaced parity test [2], whereas 
in Ref. [3] its nonlocality was investigated using the pseudospin operators jlj when only 
dissipation occurs. 

In Ref. jS] we have suggested a conditional measurement scheme on TWB leading 
to a non-Gaussian entangled mixed state, which improves fidelity of teleportation of 
coherent states. This process, termed inconclusive photon subtraction (IPS), is based 
on mixing each mode of the TWB with the vacuum in an unbalanced beam splitter and 
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then performing inconclusive photodetection on both modes, i.e. revealing the reflected 
beams without discriminat ing the number of the detected photons. IPS states have the 
following properties: they improve the teleportation fidelity for coherent states [Sj and 
show enhanced nonlocal correlations in the phase space [üj in ideal conditions, namely 
in the absence of noise. Motivated by these results and by the recent experimental 
generation of IPS states [7j, in this paper we extend the previous studies on the TWB 
and consider the nonlocality of the IPS state in the presence of noise. 

The paper is structured as follows. In Sec.|2]we address the evolution of the TWB in 
a noisy channel where both dissipation and thermal noise are present, whereas in Sec. 
we briefly review the IPS process. In Secs. EJ El and IB] we investigate the nonlocality of 
TWB and IPS by means of three different tests: displaced parity, homodyne detection, 
and pseudospin test, respectively. Finally, Sec. closes the paper with some concluding 
remarks. 



2. Dynamics of TWB in noisy channels 

The so called twin-beam state of radiation (TWB), i.e. |A)) = Vl^X* J2 k X<2 \ k ) ® \ k ) 
with À = tanhr, r being the TWB squeezing parameter. |A)) is obtained by parametric 
down-conversion of the vacuum, |A)) = exp{r(aT&' — ab)}\0), a and b being field 
operators, and it is described by the Gaussian Wigner function 

o(,w = wbsr ' () 

with 

Ao = , ° = , Bn = : , (2) 

16v/Det[cr ] 16^Det[(T ] 

where Aq = Aq(v) = cosh(2r), Bo = Bq(v) = sinh(2r) and cr is the covariance matrix 

1 / A 1 2 B o-3 \ 
^o = 7 , 3 

I2 being the 2x2 identity matrix and cr 3 = Diag(l, —1). Using a more compact form, 
Eq. (0) can also be rewritten as 

47rV D et [tr ] 



with X = (xi, yi, x 2 , y-i) T 1 cu = X\ + iyi and j3 = x 2 + iy2, and (• ■ -) T denoting the 
transposition operation. 

When the two modes of the TWB interact with a noisy environment, namely in 
the presence of dissipation and thermal noise, the evolution of the Wigner function ([!} 
is described by the following Fokker-Planck equation jHl 13 EI 

d t W t (X) = l -{d T x TX + d T x Ta 00 d x )W t {X) , (5) 
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with dx = (d Xl ,d yi ,d X2 ,d y2 ) T . The damping matrix is given by F 
whereas 





k=l 









<7 (2) 

f oo 



fc=l F ^ 1 2 



(6) 



where is the 2x2 null matrix and 

, f l + 2iV fc 



(7) 



1 + 2N k 

T k , N k denotes the damping rate and the average number of thermal photons of the 
channel k, respectively. cr^ represents the covariance matrix of the environment and, 
in turn, the asymptotic covariance matrix of the evolved TWB. Since the environment 
is itself excited in a Gaussian state, the evolution induced by © preserves the Gaussian 
form (j3J) . The covariance matrix at time t reads as follows [Tül ITT] 

(8) 



G] /2 <t ®] /2 + (i-G t )a 0O , 



where G t = ©Li e ~ Tkt la- The covariance matrix er t can be also written as 



A t (T 1 ,N 1 )l 2 B t {T 1 )<T 3 
B t (T 2 )a 3 A t (T 2 ,N 2 )l 2 



with 



A t (T k , N k ) = A e- Fkt + (1 - e~ r ^) (1 + 2N k 



B t (T k ) = B fí e- Vkt 



(9) 



(10) 



k) — I 

Let us now consider channels with the same damping rate T but different number 
of thermal photons, N\ and using the density matrix formalism, the state 

corresponding to the covariance matrix (jüj) has the following form 

Qt = 5 2 (0/ii®/i 2 5 2 t (0, (11) 
where fi k is the thermal state 



1 



1 + M k V 1 + Mi 



M k 



(12) 



a k , k — 1, 2 being the mode operators. The average number of photons are given by 
1 



Mi 
M 2 



A 2 + - 16B t - (2 - A. 



A\ - 16B t - (2 + A. 



(13) 
(14) 



A± = A u ± A %u A k , t = A t (T,N k ) and B t = B t (T). In Eq. Jü} S 2 (Ç) = exp{Ça\al - 
£,* a i a 2} denotes the two-mode squeezing operator, with parameter Ç G C 



= sinh" 1 
arg[£] =tt/2. 



-4, 



2(A\ - lQB t ) l l 2 



(15) 
(16) 
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Figure 1. Scheme of the IPS process. 



Eq. (jllj) says that the quantum state of a TWB, after propagating in a noisy channel, 
is the same of a state obtained by parametric down-conversion from a noisy background 
[TT] . Their properties, and in particular entanglement and nonlocality, can be addressed 
in an unified way using Eq. (|TT|) or, equivalently, Eqs. (jïïj) and ffTÜ|) . 

Finally, if we assume Ti = T 2 = T and Ni = N 2 = N, then the covariance matrix 
(jüj) becomes formally identical to (jïïj) and the corresponding Wigner function reads 

= exp{-2l f (|q| 2 + |/3| 2 ) + 2B t (aP + a*P*)} 
4vr 2 v ^etR 



with 



~A t= A f' N ï , g, = B -' r > , (18) 
16^/Detto-t] lGv^etK 



whereas the density matrix, mutatis mutandís, is still given by Eq. (jlljl 



3. De-Gaussification and noise 



When thermal noise and dissipation affect the propagation of an entangled state, its 
nonlocal properties are reduced and, finally, destroyed P EU E2|. Therefore it is of 
interest to look for some technique in order to preserve, at least in part, such correlations, 
or to enhance the nonlocality of the state which will face the lossy transmission 
line. Since it has been shown that the de-Gaussification of a TWB can enhance its 
entanglement in the ideal case and since non-Gaussian states can be produced using the 
current technology in this and the following Sections we will investigate whether or 
not this process can be useful also in the presence of noise. 

The de-Gaussification of a TWB can be achieved by subtracting photons from 
both modes In Ref. jHj we referred to this process as to inconclusive photon 

subtraction (IPS) and showed that the resulting state, the IPS state, can be used to 
enhance the teleportation fidelity of coherent states for a wide range of the experimental 
parameters. Moreover, in Ref. jSj, we have shown that, in the absence of any noise during 
the transmission stage, the IPS state has nonlocal correlations larger than those of the 
TWB irrespective of the IPS quantum efficiency (see also Refs. [Tïïl Hïïj). 
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First of all we briefly recali the IPS process, whose scheme is sketched in Fig. ^ 
The two modes, a and b, of the TWB are mixed with the vacuum (modes c and d, 
respectively) at two unbalanced beam splitters (BS) with equal transmissivity; the 
modes c and d are then revealed by avalanche photodetectors (APDs) with equal 
efficiency, which can only discriminate the presence of radiation from the vacuum: the 
IPS state is obtained when the two detectors jointly click. The mixing with the vacuum 
at a beam splitter with transmissivity T followed by the on/off detection with quantum 
efficiency n is equivalent to mixing with an effective transmissivity r 

t = t(T,ti) = 1-ti(1-T), (19) 

followed by an ideal (i.e. efficiency equal to 1) on/off detection. Using the Wigner 
formalism, when the input state arriving at the two beam splitters is the TWB Wo(Q) fl) 
of Eq. (JU), the state produced by the IPS process reads as follows (see Ref. for the 
details about the calculation and about the de-Gaussification map for the density matrix 
and Wigner function in the case of a TWB) 

wTVfl - 2 \ (20) 

7r 2 pii(r,r) ^ 



k=l 



where 



Pu{r,T)=y = (21) 

t^(b-fk)(b-g k )-(2B T + h k f 

is the probability of a click in both the APDs. In Eqs. ÍI20D and (|21jl we introduced 

W${a,l3) = exp{-(6 - f k )\a\ 2 - (b - g k )\(3\ 2 + (2B r + h k )(a{3 + a*(3*)} , (22) 

and defined 

C k (r, r) = - — ^ , (23) 

V5rtR[x Jfc y fc -4Bg(l-r) 3 ] 

where C 1 = 1, C 2 = C 3 = -2, C 4 = 4; x k = x k (r, r), and y k = y k (r, r) are 

xi = x 3 = y 1 = y 1 = a 

x 2 = %4 = Z/3 = Vi = a + 2 

with a = a(r, r) = 2[A (1 - r) + r], b = b(r, r) = 2[A t + (1 - r)]; finally, / fc , gf fc , and 
h k depend on r and r and are given by 

f k = Af k [x k B 2 Q + ABlil - 1 )(1 - r) + y k {l - À ) 2 ] , (24) 
9k = M k [x k {l - I ) 2 + 45 2 (1 - I )(l - r) + y fe 5 2 ] , (25) 
ftfc = A4 {(x fe + ï/ fc )Bo(l - ^o) + 2B [B 2 + (1 - 1 ) 2 ](1 - r)} , (26) 

Af k =Af k (r,T) = l ~ P ■ (27) 

a; fc2 / fc -4 J B 2 (l-r)2 

The state corresponding to Eq. (j2*ü|) is no longer a Gaussian state and its nonlocal 
properties, in ideal conditions, were studied in Ref. jü]. 
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Here we are interested in the case when the IPS process is performed on a TWB 
evolved in a noisy environment with both the channels having the same damping rate 
and thermal noise. The Wigner function of the state arriving at the beam splitters is 
now given by Eq. (fTTj) . and the output state is still described by Eq. (J2Uj) . but with the 
following substitutions 

Ao^Àt, B ^B t , o-o^vt- (28) 

We will denote with W^ P N S) (ot, 0) the Wigner function of this degraded IPS state. 

In the next Sections we will analyze the nonlocality of the IPS state in the presence 
of noise by means of Bell's inequalities. 

4. Nonlocality in the phase space 

Parity is a dichotomic variable and thus can be used to establish Bell-like inequalities 
[T7] . The displaced parity operator on two modes is defined as j2] 

n(a,/3) = D a (a)(-lf a Dt(a) ® D b ((3)(-lf b Dl(P) , (29) 

where a,/3 G C, a and b are mode operators and D a (a) = exp{<W — a*a} and D b {(3) 
are single-mode displacement operators. Since the two-mode Wigner function W(a, (3) 
can be expressed as [TT] 

W(a,(3) = ^U(a,(3), (30) 

n(a, (3) being the expectation value of IÍ(a, /3), the violation of these inequalities is also 
known as nonlocality in the phase-space. The quantity involved in such inequalities can 
be written as follows 

B DP = n(a l5 A) + n(a 2 , /3i) + U( ai ,(3 2 ) - Tl(a 2 , fa) , (31) 

which, for local theories, satisfies |Bdp| < 2. 

Following Ref. 0, one can choose a particular set of displaced parity operators, 
arriving at the following combination [ü] 

B DP (J) = ntv 7 ^, -y/J) + U(-3y/j, -Vj) + U(Vj, 3^J) - n(-3v / J, 3VJ) , (32) 

which, for the TWB, gives a maximum Í3dp = 2.32, greater than the value 2.19 obtained 
in Ref. [2]. Notice that, even in the infmite squeezing limit, the violation is never 
maximal, i.e. |£>dp| < 2v^2 [TH] . 

In Ref. jü] we studied Eq. (J32]) for both the TWB and the IPS state in an ideal 
scenario, namely in the absence of dissipation and noise; we showed that, using IPS, the 
maximum violation is achieved for r — > 1 and for vàlues of r smaller than for the TWB. 

Now, by means of the Eq. (j2Uj) and the substitutions (j2EJ), we can study how noise 
affects í3dp- The results are showed in Fig. EJ as one may expect, the overall effect 
of noise is to reduce the violation of the Bell's inequality. When dissipation alone is 
present (JV = 0), the maximum of violation is achieved using the IPS for vàlues of r 
smaller than for the TWB, as in the ideal case. On the other hand, one can see that 
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Figure 2. Plots of the Bell parameters Bbp for the TWB (top) and IPS (bottom); 

we set J = 1.6 x 10~ 3 and r = 0.9999. The dashed lines refer to the absence of noise 

(Tt = N = 0), whereas, for both the plot, the sòlid lines are Bdp with Tt = 0.01 

and, from top to bottom, N = 0,0.05,0.1, and 0.2. In the ideal case the màxima are 
b (twb) = 2 32 and fí (rps) = 2 respectively 



the presence of thermal noise mainly affects the IPS results. In fact, for Tt = 0.01 and 
N = 0.2, one has |i3j>p WB ' , | > 2 for a range of r vàlues, whereas |£> D I p S ''l falls below the 
threshold for violation. 

We conclude that, considering the displaced parity test in the presence of noise, the 
IPS is quite robust if the thermal noise is below a threshold value (depending on the 
environmental parameters) and for small vàlues of the TWB parameter r. 



5. Nonlocality and homodyne detection 

In principle there are two approaches how to test the Bell's inequalities for bipartite 
state: either one can employ some test for continuous variable systems, such as that 
described in Sec. EJ or one can convert the problem to Bell's inequalities tests on two 
qubits by mapping the two modes into two-qubit systems. In this and the following 
Section we will consider this latter case. 

The Wigner function WQ 1FS \a, (3) given in Eq. (|2T)|) is no longer positive-defmite 
and thus it can be used to test the violation of Bell's inequalities by means of homodyne 
detection, i.e. measuring the quadratures x$ and x 9 of the two IPS modes a and 
b, respectively, as proposed in Refs. [13 [Tü]. In this case, one can dichotomize the 
measured quadratures assuming as outcome +1 when x > 0, and —1 otherwise. The 
nonlocality of W^ IPS ^(a, (3) in ideal conditions has been studied in Ref. P where we also 
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Figure 3. Plots of the Bell parameter £>hd for the IPS states for two different vàlues 
of the homodyne detection efficiency: i]n = 1 (top), and t?h = 0.9 (bottom). We set 
t = 0.99. The dashed lines refer to the absence of noise (rí = N = 0), whereas, 
for both the plots, the sòlid lines are 2?hd with Tt = 0.05 and, from top to bottom, 
N = 0,0.05,0.1 and 0.2. 



discussed the effect of the homodyne detection efficiency t/h- 

Let us now we focus our attention on ( a , P), namely the state produced when 

the IPS process is applied to the TWB evolved through the noisy channel. After the 
dichotomization of the homodyne outputs, one obtains the following Bell parameter 

Í3 HD = E(&i, <pi) + ip 2 ) + E{$ 2 , <p x ) - E{$ 2 , cp 2 ) , (33) 

where and <pk are the phases of the two homodyne measurements at the modes a and 
b, respectively, and 

E{ïï h , ip k ) = / dxt h dx^ sign[a; tíh x Vk ] P{x, &h ,x ípk ) , (34) 

P(x$ h , Xtp k ) being the joint probability of obtaining the two outcomes x& h and x< Pk [16J. 
As usual, violation of Bell's inequality is achieved when |£>hd| > 2. 

In Fig. Ejwe plot £>hd for i?i = 0, ï9 2 = tt/2, <fi = —tc/A and <p 2 = tt/4: as for the 
ideal case |H1 CE] , the Bell's inequality is violated for a suitable choice of the squeezing 
parameter r. Obviously, the presence of noise reduces the violation, but we can see 
that the effect of thermal noise is not so large as in the case of the displaced parity test 
addressed in Sec. |U (see Fig. |2J. 

Notice that the high eficiències of this kind of detectors allow a loophole-free test 
of hidden variable theories |TH] , though the violations obtained are quite small. This is 
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due to the intrinsic information loss of the binning process, which is used to convert the 
continuous homodyne data in dichotomic results [2*üj . 

6. Nonlocality and pseudospin test 

Another way to map a two-mode continuous variable system into a two-qubit system 
is by means of the pseudospin test: this consists in measuring three single-mode 
Hermitian operator S k satisfying the Pauli matrix àlgebra [Sh, Sk] = 2ie hk iSi, S k = I, 
h, k, l = 1,2,3, and Shki is the totally antisymmetric tensor with £123 = +1 [HJII]. For 
the sake of clarity, we will refer to Si, S 2 and S 3 as S x , S y and S z , respectively. In this 
way one can write the following correlation function 

£(a,b) = ((a.S)(b-S)>, (35) 

where a and b are unit vectors such that 

a • S = costf a S z + sintf a [e^S. + e**>S + ) , 

b ■ S = cos d b S z + sin ï9 b (e iipb S. + e~ iipb S + ) , 1 ° (> ' 

with S± = h(S x ± S y ). In the following, without loss of generality, we set tpk = 0. 
Finally, the Bell parameter reads 

B PS = Efa, bi) + Efa, b 2 ) + Efa, bi) - Efa, b 2 ) , (37) 

corresponding to the CHSH Bell's inequality |Bps| < 2. In order to study Eq. (jïïTj) 
we should choose a specific representation of the pseudospin operators; note that, as 
pointed out in Refs. j2HI22I, the violation of Bell inequalities for continuous variable 
systems depends, besides on the orientational parameters, on the chosen representation, 
since different Sk leads to different expectation vàlues of Bps- Here we consider the 
pseudospin operators corresponding to the Wigner functions 



W x (a) = - signat*]], W z (a) = -^ (2) («), (38) 

where V denotes the Cauchy's principal value. Thanks to ()38jl one obtains 

/ u\ a q , 2sint? a sin^ fe . 
-C'TWBl.a, b) = cosv a cos itb H arctan [sinhíirj , (40) 

7T 

for the TWB, and, for the IPS, 
4 

L 



-^IPSl a , D J 



cos # a cos i?;, 2 sin $ a sin i?;, / 25 r + 
1 arctan 



nA k \ \/Aí 



(41) 



where ^ = (b — f k ){b — g k ) — (2B t + h k ) 2 , and all the other quantities have been 
defined in Sec. El 

In Fig.|3]we plot Bps for the TWB and IPS in the ideal case, namely in the absence 
of dissipation and thermal noise. For all the Figures we set i? ai = 0, $ a2 = n/2, and 
í9f,i = —^6 2 = n/A. As usual the IPS leads to better results for small vàlues of r. 
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Figure 4. Plots of the Bell parameter Bps in ideal case (Tt = N = 0): the dashed linc 
refers to the TWB, whereas the sòlid lines refer to the IPS with, from top to bottom, 
t = 0.9999, 0.99, 0.9, and 0.8. There is a threshold value for r below which IPS gives 
a higher violation than TWB. Note that there is also a region of small vàlues of r for 
which the IPS statc violates the Bell's inequality while the TWB does not. The dash 
dotted line is the maximal violation value 2^2- 

2.6 
2.5 
2.4 

B PS 2.3 
2.2 
2.1 
2 
1.9 

0.5 1 1.5 2 2.5 3 

r 





Figure 5. Plots of the Bell parameter 2?ps for Tt = 0.01: the dashed line refers 
to the TWB, whereas the sòlid lines refer to the IPS with, from top to bottom, 
t = 0.9999,0.99,0.9, and 0.8. The same comments as in Fig.Ustill hold. 



Whereas fígj as t — > oo, £>pg has a maximum and, then, falls below the 

threshold 2 as r increases. It is interesting to note that there is a region of small vàlues 
of r for which fíg WB) < 2 < B^ s \ i.e. the IPS process can increases the nonlocal 
properties of a TWB which does not violates the Bell's inequality for the pseudospin 
test, in such a way that the resulting state violates it. This fact is also present in the 
case of the displaced parity test described in Sec. El but using the pseudospin test the 
effect is enhanced. Notice that the maximum violations for the IPS occur for a range of 
vàlues r experimentally achievable. 

In Fig. El we consider the presence of the dissipation alone and vary r. We can 
see that IPS is effective also when the effective transmissivity r is not very high. We 
take into account the effect of dissipation and thermal noise in Figs. |ül and [7| we can 
conclude that IPS is quite robust with respect to this sources of noise and, moreover, 
one can think of employing IPS as a useful resource in order to reduce the effect of noise. 
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Figure 6. Plots of the Bell parameter Sps for diffcrcnt vàlues of Tt and in the absence 
of thermal noise (N = 0): the dashed lines refer to the TWB, whereas the sòlid ones 
refer to the IPS with r = 0.9999; for both the TWB and IPS wc set, from top to 
bottom, Tt = 0,0.01,0.05, and 0.1. The dash dotted line is the maximal violation 
value 2\/2. 
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Figure 7. Plots of the Bell parameter Bps for Tt = 0.01 and different vàlues N = 0: 
the dashed lines refer to the TWB, whereas the sòlid ones refer to the IPS with 
r = 0.9999; for both the TWB and IPS we set, from top to bottom, N = 0,0.01,0.1, 
and 0.2. 





7. Concluding remarks 

We have addressed three different nonlocality tests, namely, displaced parity, homodyne 
detection and pseudospin test, on TWB and IPS in the presence of noise. We have 
shown that the IPS process on TWB enhances nonlocality not only in ideal cases, but 
also when noise (dissipation and thermal noise) affects the propagation. As in the ideal 
situation, the enhancement is achieved when the TWB energy is not too high (small 
squeezing parameter r), depending on the environmental parameters. Moreover, in the 
case of the pseudospin test, we have seen that there is a region of small r for which the 
TWB itself does not violates the Bell's inequality, wheres after the IPS process it does. 

Finally, we mention that the enhanced nonlocality also in the presence of noise 
makes the IPS states useful resources for continuous variable quantum information 
processing. 
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